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Abstract 

Coordinated observations of the nearby BL Lacertae object Mrk 421 
obtained during May 1980, January 1984, and March 1984 are described. 
These observations give a time-frozen picture of the continuous spectrum 
of Mrk 421 at X-ray, ultraviolet, optical, and radio wavelengths. The 
observed spectra have been fitted to an inhomgeneous relativistic jet 
model like the one described in Hutter and Mufson (1986). In general, 
the models reproduce the data well. Many of the observed differences 
during the three epochs can be attributed to variations in the opening 
angle of the jet and in the angle that the jet makes to the line of 
sight. The jet models obtained here are compared with the homogeneous, 
spherically symmetric, synchrotron self-Compton models for this source 
computed by Makino et al . (1987). Our models are also compared with the 
relativistic jet models obtained for other active galactic nuclei. 
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I. Introduction 

In a previous paper (Mufson ejt al . 1984, Paper I), we presented 
time-frozen spectra of the BL Lacertae objects Mrk 180 and Mrk 501 which 
were obtained during a campaign of coordinated, multi frequency 
observations from the radio to the X-ray. In Paper I we fitted these 
time-frozen spectra with a variety of simple models, including the 
synchrotron-self Compton model (SSC) and a relativistic jet model 
(Hutter and Mufson 1986, Paper II). The conclusion we reached in Paper 
I was that the relativistic jet model provided the most reasonable 
description of the data sets. 

The nearby (z = .031) BL Lacertae object in Mrk 421 (B2 1101 + 38) is 
identified by its nonthermal power law continuum in the radio (Colla et 
al . 1975; Owen et al . 1978), optical (Maza, Martin, and Angel 1978; 

Muf son et al . 1980) , and X-ray (Cooke et al . 1978) ; the absence of 
emission lines in its optical spectrum (Ulrich et al. 1975); and the 
strong variability observed at all wavelengths (Miller 1975; Ricketts, 
Cooke, and Pounds 1976; Purton, Kojoian, and Dickinson 1977). The 
optical polarization (Maza, Martin, and Angel 1978) and flat radio 
spectrum (O’Dell et al . 1978) observed in this source are also common 
characteristics of BL Lacertae objects. In addition, Mrk 421 has an 
X-ray spectrum with a variable, high-energy (> 10 keV) component 
(Mushotzky et al . 1979). 

Photometric and spectroscopic observations indicate that there is an 
extended optical envelope surrounding Mrk 421 which has the properties 
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of a giant elliptical galaxy (Ulrich et al . 1975; Miller, French, and 
Hawley 1978; Mufson et al . 1980; Hickson et al . 1982)). The picture of 
Mrk 421 which emerges is that of a typical BL Lacertae object in which 
an active galactic nucleus is embedded at the core of a giant elliptical 
galaxy 

In a recent paper (Makino et al . 1987), coordinated, multifrequency 
observations of Mrk 421 during January 1984 and March 1984 were 
presented and analyzed in terms of the SSC model In this paper, we add 
a third set of coordinated, multi frequency observations of Mrk 421 from 
May 1980. We then discuss all three sets of data in terms of the 
relativistic jet model, with the goal of determining how the jet 
parameters vary with time. In addition, we compare our results to the 
SSC model fits of Makino et al . (1987), looking in particular for 
differences in these two approaches. Finally, we compare the 
characteristics of our jet models for Mrk 421, Mrk 180, and Mrk 501 with 
jet models computed by others. 

II. Observational Data 

The data sets used in the fits are from 3 epochs - May 1980, January 
1984, and March 1984. The two data sets for the 1984 observations are 
reported and discussed in Makino et al . (1987). The 1980 data set 
reported here was analyzed in the same manner as the data for Mrk 180 
and Mrk 501 reported in Paper I. The reader is referred to those papers 


for details. 
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a) Analysis of the 1980 Observations 
The 1980 data set used in the model fits is given in Table 1. A brief 
discussion of the observational results is given below. 

i ) X-Ray Regime 

The X-ray data for 1980 were obtained by the IFC and MFC aboard the 
HEAO-2 satellite (Giacconi et al . 1979). Since the IFC field in the Mrk 
421 observations has a strong, bright X-ray source within 1 arc min of 
the optical position of the active galactic nucleus, we associate the 
X-ray source with the BL Lacertae object. The absence of other strong 
X-ray sources in the IPC field argues that the MFC flux we observed also 
comes from Mrk 421. 

In Table 1 we give the details of the HEAO-2 observations. The IPC 
X-ray data reported here are the reprocessed PHASE 2 data The power 

law spectral fits to the IFC data are given in the IPC data output. In ' 

2 

Figure 1 we show contours of \ /dof in the (a, Njj) plane for these fits. 

The errors quoted in Table 1 are our 1 a estimates from this figure. 

2 

Although (\ /dof) m j n = 3.9 suggests that the fit is marginal, the 

fitting algorithm classifies it as "good”, its highest quality 

indicator. This suggests that systematic errors contribute to x • We 

consequently regard the HEAO-2 power law fit as acceptible. At 4.5 keV, 

the IPC fit gives a flux density of ~ 3 pjy which is greater thano the 

TENMA X-ray flux density observed in March 1984 and January 1984 (Makino 

90 —9 

et al . 1987). We note that Heiles (1975) measures Njj < 1.7x10 cm at 

21 cm, in good agreement with these IPC results. j 

i 

! 

i 
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The MFC data have also been fit with a power law spectrum. These 
results are also given in Table 1. In these fits the column density of 
hydrogen is unreliable since it is difficult to determine soft X-ray 
absorption with the MFC. It appears from the MPC data that the spectrum 
is steepening at higher energies. We regard the MPC spectral fits to be 
of lower overall quality than the IPC fits. 

For the purposes of our fitting algorithm as described in Paper II, we 
have broken up the TENMA spectral observations into two power laws - a 
low energy (1.5-4. 5 keV) component and a high energy (>4.5 keV) 
component . For the March 1984 observations, the X-ray spectrum appears 
to be steepening and the power laws were taken to be: log S (Jy) = -5.79 

- 1.7[log u(Hz) - 17.68] at low energies and log S (Jy) = -6.15 

- 4.9[log u(Hz) - 17.99] at high energies. The errors in the slopes are 
approximately ±0.3 and ±0.5 for the low energy and high energy power 
laws, respectively; the errors in the normalizations are approximately 
±0.02 and ±0.15 for the low energy and high energy power laws, 
respectively. For January 1984, the observations were more difficult to 
analyze. The low energy power law was taken to be of the form 

log S (Jy) = -5.47 -1.6(log u(Hz) -17.68] at low energies; the errors 
are ±0.02 and ±0.2 in the normalization and spectral index, 
respectively. The high energy data, however, are much more difficult to 
characterize by a single power law. We therefore bracketed the high 
energy data with two power laws: a "flat” power law of the form 
log S (Jy) = -6.11 - .14[log u(Hz) - 17.99] and a "steep" power law of 
the form log S (Jy) = -6.11 -2[log u(Hz) - 17.99]. 
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In Figures 2, 3, and 4, the X-ray data for the 3 epochs are plotted. 
In Figure 3, the low energy and high energy X-ray power laws are 
represented by double solid lines. The "flat" and "steep" high energy 
X-ray power laws bracket the high energy X-ray emission. 

ii) Ultraviolet Regime 

The ultraviolet observations were obtained using the IUE satellite 

(Boggess et al . 1978). Single spectral observations of Mrk 421 were 

made in the low-dispersion mode (effective resolution approximately 10 

X) with both the SWP (1150-1950 X bandpass) and LWR (1900-3200 X) 

cameras. The image numbers of the observations are given in Table 1. 

Spectral images were assembled from geometrically and photometrically 

rectified data samples in the line-by-line files of the IUE guest 

observer tapes. The long exposure, faint images of Mrk 421 were 

corrected for radiation events and fixed pattern noise as discussed in 

Paper I. The power law fit to the IUE data is given in Table 1. The 

fit listed in Table 1 is a composite of the SWP and LWR data. The 

spectral index in May 1980 appears to be flatter than it was in 1984. 

2 

The x /dof for this fit, however, is not as good as it is for the two 

1984 IUE data sets. The Integrated flux (1200 A - 3000 A) for May 1980 
-11 -2 -1 

is 5.8x10 erg cm s , which is similar to the value for January 
1984 and greater than that for March 1984 (Makino et al . 1987). 

iii) Optical Regime 

The ground based optical photometry consists of braodband BVRI 
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measurements made using the 1.5-m telescope of the University of Arizona 
Mt. Lemmon Observatory. The data were acquired using dc techniques. 
Extinction and transformation coefficients were determined from 

observations of primary UBVRI standards (Johnson et al . 1966) . The V 

-2 

band flux density of 1.7x10 Jy is signif icantly greater than any of 
the V band measurements made in January or March 1984 (Makino et al . 
1987) . 

iv) Radio Regime 

The radio data were obtained with the University of Michigan computer 
controlled 26-m radio dish. The flux densities given are means of the 
observations made that month weighted by their photometric errors. The 
errors quoted are 1 a. The data do show evidence for variability at the 
15-25% level. However, the shape of the spectrum always remains quite 
flat. Consequently, in our model fits we used the flux at the 8.0 GHz 
and flux densities at 4.8 GHz and 14.5 GHz extrapolated from the average 
spectral index observed between 1980-1982, a ~ 0. The 14.5 GHz flux 
density we used is consistent with the 14.5 GHz observed in 1984 (Malcino 
et al . 1987). 

v) Composite Spectra 

In Figures 2, 3, and 4 we have plotted the data for the 3 epochs. 

We note that although the flux in the different bandpasses has varied 
between May 1980 and January/March 1984, only the changes at optical and 
X-ray wavelengths show any evidence for correlated behavior. 
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Examination of figures 2, 3, and 4 show that both the optical and X-ray 
have declined over the three epochs of observation. 


III. Relativistic Jet Models for Mrk 421 

We have used a relativistic jet model, similar to that of Blandford 
and Konigl (1979) and Konigl (1981), to model the simultaneous radio, 
optical, ultraviolet, and X-ray observations of Mrk 421. Relativistic 
jet models like these can explain many properties of AGNs (Paper II). 
There is a great deal of observational evidence suggesting that jet-like 
features are common in AGNs (Konigl 1986). 

Details of the model and the fitting procedure can be found in Paper 
II. 

a) The Model 

In our model the relativistic motion takes the form of conical jet of 
opening half-angle ^ which is viewed at an angle 0 to its axis. It is 
assumed that the jet is free with no external pressure to confine the 
expansion in the direction perpendicular to the axis. For a free jet 
whose material is expanding with bulk velocity = Vj/c, the opening 
half-angle is given by $ 1 /Tj, where Tfj = (1-0^) ^ ^ is the Lorentz 

factor of the jet plasma. 

In the comoving frame of the expanding jet gas, the emission is 
assumed to be locally described by the homogeneous synchrotron 
self-Compton mechanism. It is further assumed that the population of 
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relativistic electrons is continuously reacclerated so that the 
distribution of electrons with respect to energy, n (nr )d*r = 
Knr^^e^^dnr , is independent of time. Continuous reacceleration is 
required since the lifetimes of the synchrotron electrons radiating at 
high energies are so short compared with the lifetime of the object. 

The wel 1-tangled magnetic field strength B and the particle number 
density K along the jet are of the form 

B = BjCr/l pc)" m , (1) 

and 

K = K^r/l pc) _n * (2) 

where r is the radial coordinate measured from the jet apex. As viewed 

by an observer on earth, the radiation is beamed in the forward 

direction by a factor Dj = nr^( 1-/3 jCosG) ^ , where nr j and /3j measure the 

bulk velocity of the jet gas. By beaming its radiation in the forward 

direction, these models avoid the problem of huge quasar luminosities 

from extremely small volumes. As discussed in Konigl (1981) and Paper 

II, each volume element in the jet contributes a local component of 

non thermal emission to the composite jet spectrum. The local jet 

spectrum can be characterized by several distinct emission regions. 

Below a frequency v , the source is self-absorbed and its flux density 

2 5 

is proportional to d ' , Below v , very little flux is contributed by 

sm 

any element of the jet to the total BL Lacertae continuum. The 
optically thin spectrum above has a spectral index a up to a 
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frequency where the spectral index abruptly increases to a^+0.5 due 

to synchrotron radiation losses. The highest energy electrons within 

the jet have a radius-independent Lorentz factor of nr ; above the 

frequency at which these electrons radiate, the spectrum drops off 

exponentially. Thus gives the frequency of the upper limit to the 

synchrotron spectrum at a particular position. We assume that the time 

evolution of y^ is slow compared with the reacceleration time; under 

these circumstances the local synchrotron spectrum remains stationary. 

The lower limit to the local synchrotron emission at a particular 

position is given by u , 

min s 

The superposition of many different regions of synchrotron emission in 
the jet produces a composite nonthermal spectrum which is flat below v < 
but which has a constantly varying spectral index for v > d^. In 
the observable radio and optical-ultraviolet wavelength regions, the 
spectrum can be approximated by two power laws: 


obs 


(») 


S ( U sM ){u/U sM ) “ sl 
S{v sll )(v/v sll )- a s2 


v - u sM C 3 ® 1 ) 

d < v < v , (3b) 

o “ “ uv v J 


where a ^ is the spectral index which applies in the radio region and 
a s2 * S t * ie s P ectra l index in the ultraviolet-optical region. The 
frequency at which the composite synchrotron spectrum begins to steepen 
from is given by the flux density at this frequency is S(i) g jj)._ 

The jet itself can be characterized by an inner radius r^ and an outer 
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radius r . As discussed in Paper II, the Inner boundary, r M , is 

determined from the frequency at which the spectral index begins to 

steepen from , the spectral index in the radio, to a^, the spectral 

index in the optical-ultraviolet. Similarly, r^ can be determined from 

the intensity of synchrotron emission at its low frequency cutoff 

U min s’ at *kis ra ditis the emission is dominated by electrons radiating 

at i) . 

min s 

Finally, the ratio of the strength of the optical nonthermal emission 
from the BL Lacertae object to the strength of the stellar optical 
emission from the giant elliptical galaxy is characterized by Ry, the 
ratio of these strengths in the V band. For sources with large values 
of Ry, the continuum is dominated by quasar emission; for sources with 
small values of Ry, the continuum has a continuum dominated by galaxy 
emission. 


b) Results of the Model Calculations 

The best-fitting models of the spectra of Mrk 421 during the three 

epochs of observation are shown in Figures 2, 3, and 4 superposed onto 

the data from the coordinated observations. The models were computed by 

varying the 9 parameters 0, $, u S(u „), a , u . ,tt , a n , and 

sM sM e min s eu sz 

Ry. The parameter a g j was held fixed at the observed value of the radio 
spectral index. In our fits we began by setting equal to the 
observed uv spectral index since the IUE spectra are presumably pure 
synchrotron emission uncontaminated by stellar emission from the galaxy. 
However, since the fits are relatively sensitive to the exact value of 
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a s2* we a ^ owe< ^ this parameter to vary within 3 ct of its observed value. 

2 

A x -statistic was then minimized to find the best fitting models. This 
2 

X -statistic was based upon the observed values of the radio flux at 

4.8, 8.0, and 14.5 GHz; the U, B, V, and R optical fluxes; and the 

normalization constants (S q ) and spectral indices (a) at UV, soft X-ray, 

and hard X-ray wavelengths. The number of degrees of freedom (dof) in 

these fits is thus 13-9 = 4. Since we have no U measurement for May 

1980, dof = 3 for this model. In these models the galaxy colors were 

fixed at the K-corrected colors of a normal giant elliptical galaxy 

(Hutter and Mufson 1981), although Hickson et al . (1982) suggest that 

the host galaxy of Mrk 421 may be redder than normal. 

In Table 2 the best-fit model parameters are given. In this table we 

have given two fits for the January 1984 data. These two fits are 

computed using the "flat" and "steep" power law representations of the 

2 

TENMA high energy X-ray spectrum, as discussed above. The x /dof for 

the January/March 1984 data sets show that these fits are acceptible. 

2 

The x /dof for May 1980 suggests that the model fit is only marginally 
adequate . 

These results can be most easily understood from the graphical 

displays of the models. In Figure 2, the superposition of the computed 

jet model onto the observations suggests that the model is adequate 

given the relatively low overall quality of the data. The poorer 
2 

X /dof, compared with the models for the other data sets, is easily 
understood from tin inspection of the figure. In particular, the 
inconsistency of the IPC and MFC data present constraints on the fitting 
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procedure which are impossible to overcome. This figure shows, however, 

that the spectrum is adequately represented by synchrotron emission from 

the radio to the X-ray. In Figure 3, both computed models are shown 

superposed onto the data. This figure demonstrates that both models are 

similar in the radio, optical, and ultraviolet bandpasses and reproduce 

the data quite well. In the X-ray, the emission from the models has 

been broken into a synchrotron component (dash-dot curve), which is 

similar for both models, and two self-Compton components (dashed 

curves), which represent fits to the "flat" and "steep" high energy 

power laws. It is clear that a self-Compton component is necessary to 

explain the high energy emission in these models. We recognize that the 

observed high energy tail could simply be an additional component of 

synchrotron emission. Nevertheless it is encouraging that the variable, 

hard spectral component (Mushotzky et al . 1979) can be explained by a 

single jet structure whose characteristics change with time. Neither of 

these models adequately account for the infrared data. But the fit 

o 

procedure we used does not include infrared data in the \ statistic. 

The model for March 1984 in Figure 4 shows how a suitable modification 

of the fit parameters could account for the infrared emission. If the 

composite spectrum were more steeply curving in the infrared and 

were at a lower frequency, then the infrared emission would simply be 

galaxy emission. In Figure 4, the superposition of the model onto the 

data shows that the model reproduces the data well in all bandpasses as 

2 

is clear from the value of \ /dof. 
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IV. Discussion 

a) Physical Parameters of Model Jets 

In Table 3, the physical parameters derived from the fits are given. 
The first two parameters characterize the degree of relativistic motion 
and Doppler beaming present in the models. The jet in Mrk 421 is 
clearly relativistic at all times. Due to small compressions/expansions 
in the jet opening angle (<()). however, the Doppler beaming factor varies. 
The Doppler beaming factor is an important characteristic of jet models 
because it explains the brightness of quasars without resorting to huge 
luminosities. Since the emission from a jet is restricted to small 
solid angles, the jet emits only a small fraction of the luminosity of a 
source radiating over 4ir. Consequently, the variations in Dj and the 
wobble in the jet with respect to our line of sight (0) easily account 
for the variability we observe in the source flux from Mrk 421. 

The four parameters m, Bj, n, and describe the magnetic field and 
the particle density as a function of radius in the jet (eq[l] and 
e< l[2]). Table 3 shows that the run of magnetic field and relaativistic 
particle density is similar for all fits. As discussed in Paper II, an 
isotropic (B^ ~ Bjj ) magnetic field should scale as B a r Table 3 
shows that the values of m are consistent with this scaling law. In 
addition, for a freely expanding conical jet, conservation of particle 
number in a free jet requires that n = 2. As can be seen in Table 3, 
conservation of particle number apparently holds in Mrk 421 during all 
observations. Finally, Konigl (1981) demonstrates that the special case 
m = 1 and n = 2 corresponds to equipartition between magnetic field and 
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particle energy. Table 3 shows that the jet in Mrk 421 is consistent 
with equipartition at all epochs. 

Table 3 further shows that the inner jet radius, r^, is extremely 

9 

small, close to the Schwarzschi Id radius for a 10 M black hole. The 

o 

implication is that the jet penetrates close to the most probable AGN 

engine. In this table we have not given a value for r , the outer 

radius of the Jet. For Mrk 421, the value of r u cannot be consistently 

determined since the low frequency synchrotron flux is dominated by 

emission from an extended, steep-spectrum halo component rather than 

emission from a jet (Weiler and Johnston 1980). The parameter 

r (u, = v ) indicates the radius at which radiation losses become 
sb su' 

important. It is the point at which the local synchrotron spectral 

index a changes to a +0.5. The frequency at which v , = v 

e e sb su 

indicates the frequency radiation losses become important. For all the 
models this occurs in the ultraviolet. The final parameter given is 
u (r M ) which is the frequency of the termination of the synchrotron 
spectrum. Any emission above this frequency arises from Compton 
emission or another component of synchrotron emission. For January 1984 
and March 1984, this occurs in the X-ray and would have been in 
principle observable. 

As discussed, one advantage of an inhomogeneous relativistic jet model 
is that it can naturally account for the variability observed in Mrk 421 
during the 3 epochs of observation. In detail, the comparison of 
Figures 2, 3, and 4 shows that there was little variation in the radio 
flux during the 3 epochs, but a decline in the optical and soft X-ray 
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f lux from May 1980 through January 1984 and March 1984. Our jet models 
reproduce this variability by changing the slope of the synchrotron 
spectrum and varying the overall level of the Compton emission. In the 
models, the weakest high frequency synchrotron spectrum occurs in 
January 1984. This occurs because the jet at this epoch has its lowest 
flow speed -Yj and highest density. Consequently, synchrotron radiation 
losses steepen the local spectral index to a ~ a g + 0.5 over a 
significant fraction of the jet’s length [small value of r ( u s ^ = w gu )]- 
In addition, the upper cutoff in the electron energy distribution (-r ) 

vll 

is lowered at low flow speeds. Both effects combine to lower the total 

synchrotron emission from the jet at high frequencies. Interestingly, 

although the low flow speed leads to a lower synchrotron flux at X-ray 

energies, the increased particle density results in a higher level of 

2 -2 

Compton emission at this epoch (Compton emission an D. ; Paper II). 

® J 

The added contribution of the Compton spectral component actually 
results in a slightly higher overall level of X-ray flux in January 1984 

over that two months later. Since the flux of the Compton emission is 

-2 

a Dj , the large variations seen in the X-ray flux (particularly in the 
hard X-ray component) are not surprising when viewed in the context of 
the jet models. 

In the radio, the lack of significant variability is to be expected 
since the emission is from extended regions of the jet plasma. For the 
January and March 1984 models, where the radio spectra particularly 
flat, the outer radius of the emitting region of the jet (r u ) is very 
large (Table 3). In these models a flat radio spectrum (® s j 0) 
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results from the combined emission from many optically thin regions of 
spectral index ~ 0.5, all of which differ in their low frequency 
turnover u . At frequencies below 1 GHz, the spectrum again begins to 
rise due to the extended halo of emission (Weiler and Johnston 1980). 

Variations in the optical region are the result of fluctuations in the 
flux of synchrotron emission beamed toward the observer. That is, 
changes in the viewing angle and opening angle of the jet add a greater 
or lesser component of nonthermal emission to the constant component of 
galaxy emission. These in turn lead to variations in R^ (Table 2) or 
variations in the color and intensity of the optical emission. The very 
blue colors of Mrk 421 (Hut ter and Mufson 1981) are the consequencies of 
large values of Ry. 

b) Comparison with SSC Models 

In Makino et al . (1987) homogeneous and spherical SSC models were fit 
to the data from January 1987 and March 1987. In many ways these SSC 
models are quite complementary to the relativistic jet models 
constructed here. As in jet models, SSC models demand a population of 
relativistic electrons with a power law energy distribution which emit a 
local power law spectrum of synchrotron radiation. In both the jet 
models and the SSC models for Mrk 421, these electrons must be 
continuously reaccelerated in situ to prevent synchrotron losses from 
removing the highest energy electrons from the assumed time-independent 
distribution of relativistic particles. In both models, the low 
frequency radio flux arises from an extended component; the optical flux 
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is due to a combination of beamed quasar emission and galaxy emission; 
and the high energy X-ray tail seen in January 1984 is the result of 
self-Compton emission. 

There are differences in the two approaches, however. In particular, 
the relativistic jet models described here fit a set of model parameters 
to the data; the SSC models only offer a range of model parameters based 
upon specific choices for the source size, the frequency of the 
synchrotron self-absorption turnover, and the flux density at the 
turnover. Often the model parameters in SSC models are relatively 
uncertain since they depend on the input parameters to very high powers. 
Further, relativistic jet models explain flux variations in a more 
natural way. Slight wobbles in the jet viewing angle and 
compressions/expansions in the jet opening angle would be expected in an 
inhomogeneous jet. Variations in source size, the explanation preferred 
in SSC models, seems somewhat more unphysical. But there is a price to 
pay for using relativistic jet models. Jet models have a great number 
of parameters (9 in our case) which must be fit to data which are often 
too poorly known to pin down the parameters adequately. 

In detail it is difficult to compare many of the parameters of the SSC 
models given by Makino et al . (1987) and the jet models described here. 
In jet models many of the source parameters vary with position in the 
jet and in SSC models many of the source parameters are uniform 
throughout the radiating volume. However, a few parameters can be 
compared directly. The Doppler factor, D ^ , which characterizes the 
degree of relativistic bulk motion in the source, is much greater in the 
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jet models than the value for the SSC models. The cooling times for the 
electrons emitting the bulk of the emission in the optical-uv is in the 
range of hours to days for both models. 

c) Comparison with Jet Models of Mrk 180 and Mrk 501 

In Paper I, relativistic jet models were fitted to time-frozen spectra 
of the nearby BL Lacertae objects Mrk 180 and Mrk 501. In Paper II 
these two jet models, along with a preliminary model of Mrk 421, were 
discussed and compared. The jet models for Mrk 421 discussed here, 
however, were constructed in a slightly different manner than those for 
the other two sources. In this paper, 9 jet parameters were varied in 
the fitting procedure, whereas in Paper I only 5 parameters [ u s jj* 

S(u „), u . , v , R.,1 were varied in the fit. Nevertheless, many of 

v sM' min s eu v J 

the individual comparisons between the sources discussed in Paper II are 
still valid here. 

Both Mrk 180 and Mrk 421 have values of D^, the Doppler beaming 
factor, which indicate that relativistic bulk motion is present in their 
jets. Mrk 421 is clearly the most relativistic; its jet is most closely 
aligned to our line of sight and its opening angle is smallest. This 
accounts for its large value of Ry and its very blue color (Hutter and 
Mufson 1981). The jet in Mrk 501, however, is pointed sufficiently far 
from our line of sight that there is no need to invoke relativistic 
beaming for this source; the colors of this source are most like a 
normal giant elliptical galaxy (Paper I). The structure of the magnetic 
fields in Mrk 421 and Mrk 501 are consistent with a locally isotropic 
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field (B^ ~ B|j) over the length of the jet. In fact this assumption was 
implicitly made in the synchrotron and Compton radiation formulae we 
used in our models. Mrk 180, however, does not appear to follow this 
scaling law. As discussed, the run of relativistic particle number with 
radius is different in Nrk 421 than it is in Mrk 180 and 501. In Mrk 
421, the most relativistic of the jets, the particle number is 
apparently conserved (n ~ 2); in Mrk 180 and Mrk 501, however, particle 
injection is apparently occurring (n < 2). This is difficult to 
understand physically and is most likely a statement about the 
imprecision with which we compute these parameters. 

As discussed in Paper II, a comparison of the energy density in the 
relativistic electrons and the energy density in the magnetic field at 
several radii shows that the magnetic energy density dominates the flow 
at small jet radii, while at r £ 1 pc the relativistic particle energy 
dominates the flow. This behavior is consistent with models for the 
formation of relativistic jets (Blandford and Payne 1982; Rees et al . 
1982; Phinney 1983) where the rotational energy of a black hole is 
extracted and transferred to the relativistc particles by the magnetic 
field. Near the black hole engine, where the transfer is taking place, 
the magnetic field would carry the most significant fraction of the 
energy. Farther out, the energy has been transferred to the particles 
and they would carry the bulk of the energy. Since this type of 
hydromagnetic jet acceleration is thought to occur when the accretion 
rate onto the central black hole is low (Blandford, Begelman, and Rees 
1984), it is most likely to be operating in low power AGN’s like Mrk 
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180, Mrk 421, and Mrk 501. It Is interesting to note that in the SSC 
models of Makino et al . (1987) the magnetic energy density dominates the 
relativistic particle energy. 

d) Further Comparison of Relativistic Jet Models 

Urry and Mushotzky (1982) and Worrall et al . (1984a, b.c; 1986) have 
constructed relativistic jet models for many sources. In all these 
models, the authors conclude that relativistic jets provide a 
satisfactory fit to the data. Values of magnetic field strength, 
particle ennsities, and Doppler factors cover ranges similar to those 
peported in Papers I and II. 

V. Conclusions 

Coordinated radio, optical, ultraviolet, and X-ray observations of teh 
BL Lac object Mrk 421 have been presented. These observations give a 
time-frozen picture of the continuous spectrum of Mrk 421 at each of 
three epochs. The observed spectra have been fitted with an 
inhomogeneous relativistic jet model which generally reproduces the data 
well. Many of the observed differences between the spectra at the 
different epochs can be attributed to slight variations in the jet 
direction and speed which can account for the integrated intensity 
variations through changes in the shaped of the composite synchrotraon 
spectrum of the jet and (in the X-ray) through the level of the 
self_Compton emission. This is the same conclusion previously found 
from the application of this model to the BL Lacertae objects Mrk 180 



23 


and Mrk 501. Further, the differing degrees of variability in the 
individual bandpasses also arise naturally in this model. Models 
similar to the one presented here have beenb shown by other authors to 
give similarly consistent discriptions of the emitting regions in other 
blazars. The model’s principal disadvantage lies in the large number of 
parameters which must be specified on the basis of often limited data. 
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Table 1 

1980 Observations of Mrk 421 


a) X-Ray Observations: 


UT Date 

Count 

Range 

Spectral 

log N h 

iog S o 


Rate 


Index 

(cm -2 ) 

(@ 2 keV) 


(cts/s) 

(keV) 

(«) 


(Jy) 


IPC: 


1980 

12-17 May 


1.064.05 


.2-4.0 -1.24.2 (9.541 .8)xl0 19 (7.64.6)xl0 6 3.9 1 


MPC: 


2.454.09 1.3-8. 2 -2.04.5 


(843)xl0 6 0.3 


b) Ultraviolet Observations: 

UT Date Camera Exposure 

Spectral 

log S Q (mjy) 

X^/dof 

Time 

Index 

(©log v = 15.2) 


(sec) 

(a) 




1980 
17 May 


SWP (9025) 12,600 

LWR (7779) 5,400 


-0.634.2 


-2.44.02 1.5 
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c) Optical Observations : 


UT Date 

Aperture 

(arcsec) 

B 

(mag) 

V 

(mag) 

R 

(mag) 

I 

(mag) 

1980 May 12-18 

16.8 

13.96L0.2 

13.44i.02 

12.82i.02 

12.21i.02 



d) Radio Observations: 


UT Date 

s 4 . 8 (Jy> 

S 8.0 (Jy) 

s 145 (Jy) 

1980 May 

— 

0.62±.03 

— 

1980 Nov 

0.74±.06 

0.66±.03 

0 . 58± . 05 

1981 Nov 

0.73±.05 

0.76i.03 

0 . 83i . 04 

1981 Dec 

0 . 76± . 03 

0.81i.02 

0.76±.02 

1982 Mar 

0 . 83± . 04 

0.79±.05 

0.78±.02 

<S> 

0.77±. 05 

0.74i.03 

0.76±.03 


HEAO-2 processing algorithm specifies fit as "good”. 
All observations with Mt. Lemmon 1.5 m telescope. 
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Table 2 


Relativistic Jet Model Fits to Mrk 421 for 3 Epochs 



Mav 1980 

Jan 1984/flat 

Jan 1984/s teep 

March 1984 

0 (rad) 

.04 

.15 

.09 

.06 

(j) (rad) 

.03 

.12 

.08 

.06 

" S M < Hz > 

4.7xl0 12 

1.4xl0 13 

13 

2. lxlO 1,3 

3.5X10 1 ' 

s(« sM ) (Jy) 

.45 

.19 

.16 

.31 

a 

e 

.44 

.40 

.48 

.47 

v . (Hz) 

min s v ' 

1.3xl0 9 

5.8xl0 5 

8.8xl0 4 

1 . lxlO 4 

T r 

eu 

1.7xl0 5 

4.2xl0 4 

3.5xl0 4 

5.7xl0 4 

2 

a sl 

0.3 

.1 

.1 

.1 

a s2 

.93 

1.0 

1.1 

1.0 

«v 

2.4 

4.2 

4.1 

1.6 

X 2 /dof 

7/3 

6.5/4 

6.1/4 

3.5/4 


notation: after Hutter, D.J., and Mufson, S.L. 1986, An . J . . 301, 50. 
held fixed in these fits at the observed radio spectral index. 
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Table 3 



Derived Physical Parameters 

of Jets 



Mav 1980 

.Tan 1984/flat 

Jan 1984/s teen 

March 1984 

r J 

38 

8 

13 

18 

D J 

22 

7 

10 

16 

m 

.85 

.79 

CD 

.96 

Bj (G) 

.22 

.06 

.04 

.06 

n 

1.9 

2.1 

2.1 

1.9 

K x (cm -3 ) 

.16 

.33 

.47 

.26 

r M < pc > 

9.7xl0~ 4 

7 . 4xl0~ 5 

7.2xl0~ 5 

l.xlO -3 

r ^ U sb “ U su^ 

(Pc) 





13 

2.1 

2.3 

2.4 

v , = v (Hz) 

sb su v ' 

6.2xl0 15 

l.QxlO 15 

1.5xl0 15 

5.4xl0 15 

“su (r M> < Hz > 

2. lxlO 20 

3.5xl0 18 

5.4xl0 18 

9.5xl0 18 


1 


notation: after Hutter, D.J., and Mufson, S.L. 1986, Ap..T. . 301, 50. 
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Figure 


Figure 


Figure 


Figure Captions 


2 

1: Contours of \ /dof in the (spectral index, hydrogen column 

density) -plane for power law fits to the IPC X-ray data. The 
2 

X /dof data were generated by the HEA0-2 fitting algorithm. 

2: The multi frequency spectrum of Mrk 421 from the observations 

of May 1980. The data plotted are given in Table 1. The 
radio and optical data are shown as diamonds; the power law 
fits to the IUE, IFC, and MFC data are shown as solid lines 
with error bars indicating 1 a errors in the flux density at 
the ends of the bandpass. Superposed onto the data as a 
dotted line is the best fitting relativistic jet model. The 
component of emission due to self-Compton emission is shown 
as a dashed line. The fit parameters used to describe the 
jet are given in Table 2. 

3: The multifrequency spectrum of Mrk 421 from the observations 

of January 1984. The data plotted are taken from Makino et 
al . (1987). The radio and optical data used in the fits are 
shown as diamonds. Infrared data and addtional radio data 
not used in the fits are shown as X’s. The power law fit to 
the IUE data is shown as a solid line; 1 a errors at the end 
of the bandpass are indicated. The TENMA observations have 
been broken into a single low energy and two high energy 
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("flat” and "steep") power laws, as described in the text. 

The low energy X-ray power law and the "flat" and "steep" 
high energy X-ray spectra are represented by a double solid 
line. The approximate 1 a errors in the low energy component 
are indicated. Superposed onto the data is the best fitting 
relativistic jet models for both the "flat" and the "steep" 
spectrum. At the high energy end, the model flux density has 
been broken into synchrotron (dash-dots) and self-Compton 
(dashes) components. The fit parameters used to describe 
both the "flat" and "steep" jets are given in Table 2. 

Figure 4*. The multifrequency spectrum of Mrk 421 from the observations 
f March 1984. The multifrequency spectrum of Mrk 421 from 
the observations March 1984. The data plotted are taken from 
Makino et al . (1987). The radio and optical data used in the 
fits are shown as diamonds. Infrared data and addtional 
radio data not used in the fits are shown as X’s. The power 
law fits to the IUE and TENMA data are shown as solid lines; 

1 a errors at the end of the bandpass are indicated. 
Superposed onto the data as a dotted line is the best fitting 
relativistic jet model. The component of emission due to 
self-Compton emission is shown as a dashed line. The fit 
parameters used to describe the jet are given in Table 2. 


